qcno_dB = 45;
N=10000;
stdn = 8;
stdn_IQ = sqrt(stdn^2* N/2);

ro1=0.75;
ro2 = 0.5;
Dp=stdn_IQ^2; Dpe=ro1*stdn_IQ^2; Del=ro2*stdn_IQ^2;

L=chol ([Dp Dpe Dpe;
        Dpe Dp Del;
        Dpe Del Dp])';
    Nj=1000000;
    nIp= nan (1,Nj);
    nIe = nan(1,Nj);
    nIl = nan(1,Nj);
    for j =1:Nj
        nI = L*randn(3,1);
        nIp(j)=nI(1);
        nIe(j)=nI(2);
        nIl(j)=nI(3);
    end
    fprintf('Corrcoeff nIp nIp=%f\n', mean(nIp.*nIe/std(nIp)/std(nIe)));
    fprintf('Corrcoeff nIl nIp=%f\n', mean(nIl.*nIe/std(nIl)/std(nIe)));